Universal scaling, beta function, and metal-insulator transitions 
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We demonstrate a universal scaling form of longitudinal resistance in the quantum critical region 
of metal-insulator transitions, based on numerical results of three-dimensional Anderson transitions 
(with and without magnetic field), two-dimensional quantum Hall plateau to insulator transition, 
as well as experimental data of the recently discovered two-dimensional metal-insulator transition. 
The associated reflection symmetry and a peculiar logarithmic form of the beta function exist over 
a wide range in which the resistance can change by more than one order of magnitude. Interesting 
implications for the two-dimensional metal-insulator transition are discussed. 
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The scaling theory El predicted that the noninter- 
acting electrons are always localized in two-dimensional 
(2D) disordered systems. Recently, a new scaling argu- 
ment H was put forward in order to accommodate the 
newly-found 2D metal-insulator transition (MIT) || in 
zero magnetic field (B = 0), where the Coulomb interac- 
tion presumably becomes very important B. Although 
the microscopic mechanism remains unclear 0-o], with- 
out violating any general scaling principles the authors 
assumed the following leading behavior of the "beta func- 
tion" (3(g) — d[ln(g)]/d[ln(L)] for large conductance g at 
a finite length scale L: 



0(g) = (d - 2) + A/g a + 



(1) 



in which A becomes positive in the aforementioned B = 
MIT systems y,0,|), leading to a metallic phase (0 > 0) 
at the dimensionality d = 2. 

Since 0(g) < at small g (localized region), the beta 
function is then no longer a monotonic function and has 
to change sign at some finite g = g c , which corresponds 
to a quantum critical point. Experimental measurements 
have indicated || an exponential form for the conduc- 
tance with a peculiar reflection symmetry relating the 
conductance and the resistance on both sides of MIT, 
which implies 12] the following logarithmic form of the 
beta function in the quantum critical region (QCR): 



(3(g) = -ln(g/g c ). 



(2) 



In particular, v here is the correlation length exponent, 
and (0) holds at a wide range (1/4 < g/g c < 4) far 
beyond a simple small variable expansion around g = g c . 
The logarithmic form of the beta function (0) looks 
quite remarkable. Recall that in strong localization limit 
one may find, exactly, 0(g) — ln(g) — constant. But the 
inverse exponent 1/v does not show up in front of ln(g) 
as in (0) and the corresponding behavior of g should be 
quite different. So far there still lacks a good theoretical 
understanding of (0) from a microscopic model. Nev- 
ertheless, one may ask an equally important question: 
whether (0) is a property of the beta function unique 



for the B — 2D MIT system or it actually represents a 
generic scaling behavior of quantum phase transitions in- 
cluding other MIT systems with different symmetry and 
dimensionality. Unfortunately, so far there is no direct 
anwser to this question as how the scaling function of con- 
ductance behaves and what is the form of the beta func- 
tion in the QCR of various MIT systems are not known, 
although a lot of efforts have been focused on the critical 
conductance and exponent within each universality class 

In this paper, we present direct numerical evidence 
showing that the beta function (g) in fact holds for the 
following systems as well: three-dimensional (3D) Ander- 
son transitions with and without magnetic field (repre- 
senting orthogonal and unitary classes, respectively); the 
2D electrons in strong magnetic filed, i.e. the quantum 
Hall effect (QHE) system. Strikingly, v0(g) = ln(g/g c ) is 
found to be a universal function in the QCR where g/g c 
may change up to two orders of magnitude. Correspond- 
ingly the resistance is of an exponential form p xx oc e~ s 
with s = ±(co£/£) 1 / iy which also implies a reflection sym- 
metry in the same region (here £ is correlation length, 
and Co is a non- universal dimensionless constant ~ 0(1)). 
Thus (2) may well represent a "super" universality prop- 
erty associated with general quantum phase transitions. 
Furthermore, deep into the metallic region, the beta func- 
tion shows distinct behavior depending on how the resis- 
tance p xx deviates from the exponential form: 3D MITs 
and the B — 2D MIT experimental data seem to be- 
long to the same group where d[ln(l/ p xx )]/ds is a mono- 
tonically decreasing function of scaling variable s in the 
whole scaling region; on the other hand, the QHE sys- 
tem falls into a different group where d[ln(l/p xx )]/ds be- 
comes a monotonically increasing function of s. Interest- 
ingly, the experimental data of superconductor-insulator 
transition p_2| also fall into the second group, in accord 
with the speculation [|j| that the MIT in the QHE and 
superconductor-insulator transition may belong to simi- 
lar universality class. 

We consider disordered electron systems described by 
the Anderson Hamiltonian Ml: 



H 



<ij> 



H.c. + 2J uiicfci 



where the hopping integral is taken as the unit, and cj 
is a fermionic creation operator with < ij > referring 
to the nearest neighboring sites. A uniform magnetic 
flux per plaquette (along z direction) can be imposed 



by requiring </> = J^ : 



2ir/M, where the summation 



runs over four links around a plaquette in the x — y plane. 
Wi is a random potential uniformly distributed between 
(-W/2,W/2). 

We first study the 3D electron system without mag- 
netic field (dij — 0) which belongs to the orthogonal 
class. The longitudinal conductance G xx is calculated 
using Landauer formula Q. By changing the disorder 
strength W, a metal-insulator transition is found at a 
critical disorder strength W c = 16.5 [[10| at the Fermi en- 
ergy Ef = 0, with a critical conductance G c — 0.37 (in 
units of e 2 /h) and correlation exponent v — 1.6. All the 
data at different sample sizes (L = 8, 10, 12, 14 and 16) 
can be then collapsed onto two branches as a function of 
L/£ as shown in Fig. 1(a) (o curve). Note that the re- 
sistance data are plotted in the figure. (More than 2000 
configurations are taken in the average for L = 16, and 
more for smaller L's.) A 3D MIT is similarly obtained 
in the presence of strong magnetic field (unitary class). 
We have chosen two different flux strengths tfi — 2tt/M: 
M = 5 at sample sizes L = 10 and 15; and M = 4 
at sample sizes L = 8, 12, and 16, respectively. All the 
longitudinal resistance data with different </>'s and L's 
again can be scaled onto two branches (+ curve in Fig. 
1(a)). At the critical point, W c = 18.3, v = 1.43 0] and 
G c = 0.294 at E f = 0. 

Since the universality of the MIT in the unitary class is 
distinct from the one of the orthogonal class as expected 
in the scaling theory pJlCfl, two scaling curves shown in 
Fig. 1(a) are generally different from each other. How- 
ever, if we re-plot the data in terms of the scaling variable 



a = ± (co VO 



l/V 



(3) 



where the sign +(— ) corresponds to the metallic (insu- 
lating) branch, two curves of longitudinal resistance in 
the QCR can be precisely scaled together as shown in 
Fig. 1(b). Here the dimensionless constant cq has the 
non-universal values 2.27 and 1.82 for orthogonal and 
unitary class, respectively. As shown in the insert of Fig. 
1(b), the resistance in the transition region well follows 
a simple exponential form 



Pxx/Pc = exp(-s) 



(4) 



over a rather broad region: —2 < s < 2 or 1/8 < 
Pxx/Pc < 8. Because of such a wide range of s (in- 
stead of a small parameter expansion), the exponen- 
tial behavior appears very robust. In the same re- 
gion, one always finds the so-called reflection symmetry: 




L/§ 




FIG. 1. (a) Finite-size resistance ratio p xx /p c as a scaling 
function of L/£ for 3D Anderson transitions in the absence (o) 
and presence (+) of magnetic field, (b) p xx /pc as a function 
of the scaling variable s — ±(coL/£) • The insert: the 
enlarged quantum critical region where p xx /pc — e~ s . 



Pxx{s)/ p c — Pel Pxx(-s) between the metallic and insu- 
lating branches. 

Now let us consider a qualitatively different MIT in 
the QHE system where 2D electron gas is subjected to 
a strong magnetic field. By tuning the Fermi energy (or 
the density of electrons) near the lowest Landau Level 
(LL), an insulator to metal transition can be induced 
which is characterized by a one-parameter scaling the- 
ory pd| ] with an exponent v = 7/3 and p c — 1 (in 
units of h/e 2 ) |L6|. The Hall conductance here is cal- 
culated using Kubo formula. By going to large sample 
sizes (L = 24,32,48,56,64), we were able to obtain the 
scaling behavior of p xx for the QHE systems. In Fig. 
2, p xx is plotted as a function of scaling variable s de- 
fined in (pi). Again the resistance exhibits the same ex- 



ponential dependence p xx = exp(-s) (solid line) in the 
critical region covering a similar wide range of resistance 
(1/5 < p xx /p c < 5) as in the 3D MITs. In Fig. 2, 
two different disorder strengths, W — 1 and W = 4, are 
considered which represent weak and strong LL coupling 
limit, respectively. The corresponding scaling functions 
start to deviate from the exponential form beyond the 
critical region and simultaneously become IF-dependent 
in the insulating region. 

We would like to point out an interesting reflection 
symmetry for the W — 1 case: as shown in the insert 
of Fig. 2, p xx (s) on the insulating side and l/p xx (— s) 
on the metallic side perfectly coincide with each other 
over the whole scaling region and covering a resistance 
range 1/100 < p xx < 100 which is way beyond the 
critical region. Our interpretation is that at weak dis- 
order (W = 1), the particle-hole symmetry is still ap- 
proximately maintained near the lowest LL such that 
the Hamiltonian is self-dual []13|,[17| in the Chern-Simon 
boson language, which then leads to the wide range of 
the reflection symmetry. By contrast, when disorder is 
strong and all the LLs are coupled together without the 
particle-hole symmetry, the reflection symmetry only ex- 
ists around the QCR where the exponential behavior (Q) 
is followed. 

As demonstrated by the above numerical calculations, 
the scaling function of longitudinal resistance shows a 
universal exponential behavior over a wide range in 3D 
and QHE MITs. In Fig. 3, these data are plotted to- 
gether with the experimental data obtained in the B = 
2D MIT in the Si sample |]. Note that the experi- 
mental data were measured at finite temperature so the 




length scale L should be replaced by the dephasing length 
Li n oc T~ x > z . (Here z = 1 is the dynamical exponent). 
The correlation length £ oc I/To oc |<5n| - " in the transi- 
tion region H (S n is the electron density measured from 
the critical point). So scaling variable s in this case be- 
comes ±(coT'o/T) 1 / iy with v = 1.6 and cq is a dimension- 
less constant. In Fig. 3, a universal scaling function of 
the longitudinal resistance is clearly shown for all these 
systems in the QCR (with —1.5 < s < 1.5) despite their 
different symmetry classes, dimensionalities, and micro- 
scopic mechanisms of the MIT. The corresponding beta 




FIG. 2. pxx/pc as a function of the scaling variable s for the 
QHE systems both at weak (o, W — 1) and strong (+, W — 4) 
disorder strengths. The insert: the reflection symmetry of p xx 
at W = 1. 



FIG. 3. Experimental data |J of B = MIT (•) are plotted 
together with the resistances calculated for the systems of 3D 
Anderson transition (o and +) and the QHE ( A). The corre- 
sponding P functions are shown in the insert which precisely 
follow the logarithmic form (0) in the QCR. 

functions for those systems are shown in the insert of 
Fig. 3 if we define g = l/p xx , where a straight dashed 
line represents the logarithmic form of (]2|) which can be 
obtained straightforwardly from ([|) . Note that the beta 
function is multiplied by the critical exponent v in the 
insert such that the resulting function becomes universal 
in the QCR. 

Furthermore, we would like to comment on an inter- 
esting trend in the metallic region for those systems. 
In weak disorder limit of 3D MITs, the resistance ap- 
proaches to zero in power law: p xx ~ (£/£) = cqS~ v . 
The curve for B = 2D MIT system follows very closely 
to the ones of 3D MITs on the same side of the solid 
line in Fig. 3 as it drops to zero slower than exp{— s). 
In contrast, in the QHE system, the resistance devi- 
ates the solid line on the opposite side which means 
it approaches to zero even quicker than in the QCR. 
This behavior can be easily seen in its asymptotic form: 
Pxx ~ &xx ~ exp(— s u /cq) at large s limit (since in 
the QHE plateau region electrons are also localized such 



that at large L/£ limit <j xx ~ exp(—L/^), a xy — 1 and 
Aez ~ 0xx)- One may then define a generalized dimen- 
sionless function as follows: 



/3i = d[to(p)]/ds = i//3/s 



(5) 



As shown in Fig. 4, all the data fall onto the straight 
line with (3\ = 1 in the QCR. In the metallic phase the 
distinctive large-s behavior of 0i separates the metal- 
lic regime into two regions. METAL denotes the region 
where j3\ scales to zero which is followed by the 3D MITs 
as well as the experimental data of the B — 2D MIT 
system. On the other hand, ftx for the MIT in the QHE 
system diverges to infinity at s — > oo, which is denoted as 
METAL(B) region known as the "Bose" metal following 
the theoretical description |13|] . /3i for disorder-tuned 
superconductor-insulator transition jlj] in the metallic 
region is also plotted in Fig. 4 (* curve) which indeed 
shows a quick increase like in the QHE system. (Here 
the scaling variable s is of the same form used in plotting 
the experimental data of B=0 2D MIT.) According to 
Ref. (p|, these two systems should belong to the same 
category as classified by the "Bose" metal here. 
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FIG. 4. The function /3i defined in (|5j) as a function of the 
scaling variable s for those systems shown in Fig. 3. In addi- 
tion, the data for a disorder-tuned superconductor-insulator 
transition [R2| in the metallic regime are also shown for com- 
parison (*). 

We conclude by making several comments on the na- 
ture of the B = 2D MIT systems based on the present 
work. First, no matter what the microscopic mecha- 
nism is, such a 2D MIT seems to belong to a quantum 
phase transition instead of a classical phase transition (or 
crossing over): The experimental data of the resistance 
precisely coincides with those of other known MITs in 
the QCR, plotted as a function of the scaling variable 
(Lin/ d) 1 /", which covers a range of the resistance by 
more than one order of magnitude. Second, the reflection 
symmetry of resistance [p| is the natural consequence of 
the universal resistance scaling in the QCR. Finally, the 



metallic phase behaves more like a normal metal than a 
superconductor as revealed by the classification based on 
the j3\ function shown in Fig. 4. 
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